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Abstract. Recently, the author and collaborators have developed a system- 
atic program for proving the existence of homoclinic orbits in partial differen- 
tial equations. Two typical forms of homoclinic orbits thus obtained are: (1). 
, transversal homoclinic orbits, (2). Silnikov homoclinic orbits. Around the 

■ transversal homoclinic orbits in infinite dimensional autonomous systems, the 
• ' author was able to prove the existence of chaos through a shadowing lemma. 

' Around the Silnikov homoclinic orbits, the author was able to prove the exis- 

, fence of chaos through a horseshoe construction. 

Very recently, there has been a breakthrough by the author in finding 
Lax pairs for Euler equations of incompressible inviscid fluids. Further results 
have been obtained by the author and collaborators. 

> 

o 

■ 1. Introduction 

<N 

fS| ' Unlike chaos in finite dimensional systems, the area of chaos in partial differen- 

, tial equations has been a long standing open field. During the past ten years, there 

' have been some remarkable developments in the area of chaos in partial differential 

equations by the author and collaborators. Our main results can be briefly de- 
' scribed as follows: A systematic program is developed for proving the existence of 

, homoclinic orbits for perturbed soliton equations. The program involves machiner- 

^ ' ies from integrable theory, dynamical systems, and partial differential equations. 

• • , The types of homoclinic orbits thus obtained can be either transversal or of Sil- 

^ ^ ' nikov. With regard to transversal homoclinic orbits, shadowing lemmas are utilized 

, or developed to prove the existence of chaos. With regard to Silnikov homoclinic 

' orbits, Smale horseshoes are constructed to prove the existence of chaos. The main 

machineries from integrable theory are the isospectral theory and Darboux transfor- 
mations. The main machineries from dynamical systems are persistence of invariant 
manifolds, Fenichel fibers, Melnikov measurement, and other measurements. The 
main machineries from partial differential equations are local well-posedness and 
regularity. 

Another exciting development happened to Euler equations of incompressible 
inviscid fluids. Lax pairs for both 2D and 3D Euler equations was found by the 
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author and Steve Childress. A Darboux transformation for the Lax pair of 2D 
Euler equation was also found by Artyom Yurov and the author. The philosophi- 
cal significance of the existence of Lax pairs for Euler equations is fundamentally 
important. If one defines integrability of an equation by the existence of a Lax 
pair, then both 2D and 3D Euler equations are integrable. More importantly, both 
2D and 3D Navier-Stokes equations at high Reynolds numbers are near integrable 
systems. Such a point of view changes our old ideology on Euler and Navier-Stokes 
equations. 

The most important application of the theory on chaos in partial differential 
equations in theoretical physics will be on the study of turbulence. Our approach of 
studying turbulence is different from many other studies in which one starts with 
Stokes equation to prove results on Navier-Stokes equations for small Reynolds 
number. In our studies, we start with Euler equations and view Navier-Stokes 
equations at large Reynolds number as singular perturbations of Euler equations. 
For this goal, we chose the 2D Navier-Stokes equations under periodic boundary 
conditions to begin a dynamical system study on 2D turbulence. We started with a 
simple fixed point of 2D Euler equation, and studied the linearized 2D Euler equa- 
tion at the fixed point. A complete spectral theorem is obtained. In particular, we 
found unstable eigenvalues. Then, naturally we are interested in the corresponding 
unstable manifold. Such an unstable manifold forms the prototype of the attractor 
of Navier-Stokes equations at large Reynolds number. And long term turbulence 
lives around the attractor. Finding the unstable manifold is not successful yet. 
Partial success is accomplished in finding the corresponding unstable manifold for 
some Galerkin truncation. Here the unstable manifold is the 2D surface of an 
ellipsoid. Stable and unstable manifolds together form a lip-shape configuration. 
Whether or not the stable and unstable manifolds for 2D Euler equation have the 
same topology as above is not known. In fact, the existence of stable and unstable 
manifolds for 2D Euler equation has not been proved. 

2. Homoclinic Orbits 

In terms of proving the existence of a homoclinic orbit, the most common tool 
is the so-called Melnikov integral method |61| |4]. This method was subsequently 
developed by Holmes and Marsden |23) . and most recently by Wiggins 79 . For 
partial differential equations, this method was mainly developed by Li et al. |48j 
|40| jSO] |46) . 

There are two derivations for the Melnikov integrals. One is the so-called 
geometric argument |23| [75 50 48 40 . The other is the so-called Liapunov- 
Schmitt argument |12| The Liapunov-Schmitt argument is a fixed-point type 

argument which directly leads to the existence of a homoclinic orbit. The condition 
for the existence of a fixed point is the Melnikov integral. The geometric argument 
is a signed distance argument which applies to more general situations than the 
Liapunov-Schmitt argument. It turns out that the geometric argument is a much 
more powerful machinary than the Liapunov-Schmitt argument. In particular, the 
geometric argument can handle geometric singular perturbation problems. I shall 
also mention an interesting derivation in 4 . 

In establishing the existence of homoclinic orbits in high dimensions, one often 
needs other tools besides the Melnikov analysis. For example, when studying or- 
bits homoclinic to fixed points created through resonances in {n > 4)-dimensional 
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Figure 1. The graph of 

near-integrable systems, one often needs tools like Fenichel fibers, as presented in 
previous chapter, to set up geometric measurements for locating such homoclinic 
orbits. Such homoclinic orbits often have a geometric singular perturbation nature. 
In such cases, the Liapunov-Schmitt argument can not be applied. For such works 
on finite dimensional systems, see for example |32j |33| | 5U| . For such works on 
infinite dimensional systems, see for example |48| j40| . 

2.1. Silnikov Homoclinic Orbits in Nonlinear Schrodinger (NLS) Equa- 
tion Under Regular Perturbations . Consider the regularly perturbed nonlin- 
ear Schrodinger (NLS) equation 48 , 

(2.1) tqt = + 2[|g|2 - uj^]q + ie[dlq -aq + P], 

where q = q{t, x) is a complex-valued function of the two real variables t and x, 
t represents time, and x represents space. q{t, x) is subject to periodic boundary 
condition of period 27r, and even constraint, i.e., 

q{t, X + 2n) — q{t, x) , q{t, ~x) — q{t, x) . 

a; is a positive constant, a > and (3 > are constants, is a bounded Fourier 
multiplier, 

N 

9lq = " X! ^^^^qk cos kx , 
fc=l 

^fe = 1 when k < N, = 8fc^^ when k > N, for some fixed large N, and e > is 
the perturbation parameter. The following theorem was proved in |48j . 

Theorem 2.1. There exists a Eq > 0, such that for any e G (0, £0)7 there exists 
a codimension 1 surface in the external parameter space (a, (3, uj) G M+ x R+ x M+ 
where to S (^, 1), and au) < (3. For any {a,/3,uj) on the codimension 1 surface, the 
regularly perturbed nonlinear Schrodinger equation \2.1}) possesses a symmetric pair 
of Silnikov homoclinic orbits asymptotic to a saddle Qe ■ The codimension 1 surface 
has the approximate representation given by a = 1/k{uj), where k(w) is plotted in 
Figure Q 

To prove the theorem, one starts from the invariant plane 

n^{q\ a,<z-o}. 

On n, there is a saddle = vTe** to which the symmetric pair of Silnikov 
homoclinic orbits will be asymptotic to, where 

(2.2) I = - el-^/p^^^ + ■ ■ ■ , cose^^, 0e(O,|). 
Its eigenvalues are 

(2.3) A± = -e[a + ^^n'] ±2^{^+l.^- /)(3/ - - y) , 

where n = 0, 1, 2, • • • , w S (i, 1), ^„ = 1 when n < N, = 8n~^ when n > N, 
for some fixed large N, and / is given in (|2.2I) . The crucial points to notice are: 
(1). only and have positive real parts, Re{A(j"} < Re{A^}; (2). all the other 
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eigenvalues have negative real parts among which the absolute value of Re{Aj} = 
Re{A^} is the smallest; (3). |Re{Aj}| < Rc{Aq }. Actually, items (2) and (3) are 
the main characteristics of Silnikov homoclinic orbits. 

The unstable manifold W"{Qe) of Qs has a fiber representation _48_. The 
Melnikov measurement measures the signed distance between W^IQ^) and a center- 
stable manifold Wjf* |48| . By virtue of a Fenichel Fiber Theorem 48 , one can show 
that, to the leading order in e, the signed distance is given by the Melnikov integral 

M ^ / [dqF^{qo{t)){dlqo{t) - aq^{t) + (3) 



+dqFi{qo{t)){diqo{t) - aqo{t) + P)]dxdt, 

where qoit) is a homoclinic orbit for the integrable NLS, which can be generated 
through Darboux transformations '49' ; and dqFi and dqF^ are Melnikov vectors 
which can be generated through Floquet discriminant in the isospectral theory of 
the integrable NLS |49) . The zero of the signed distance implies the existence of 
an orbit in n W^^. The stable manifold W''{Q^) of Qe is a codimension 

1 submanifold in W^'^. To locate a homoclinic orbit, one needs to set up a second 
measurement measuring the signed distance between the orbit in W'^{Qs) H W^^ 
and ^{Qe) inside W^'^. To set up this signed distance, first one can rather easily 
track the (perturbed) orbit by an unperturbed orbit to an 0{e) neighborhood of 
n, then one needs to prove the size of W^{Qe) to be 0{e'^) {v < 1) with normal 
form transform. To the leading order in e, the zero of the second signed distance 
is given by 

auj{Aj) 
2 sm 

where A7 = —4190 and do is a phase shift of qo{t). To the leading order in e, the 
common zero of the two second signed distances satisfies a = 1/k(u;), where k{uj) 
is plotted in Figure ^ Then the claim of the theorem is proved by virtue of the 
implicit function theorem. For rigorous details, see 



2.2. Silnikov Homoclinic Orbits in NLS Under Singular Perturba- 
tions . Consider the singularly perturbed nonlinear Schrodingcr equation |4U) . 

(2.4) iqt = q^^ + 2[\q\^ - uj^]q + ie[q.^^ - aq + (3] , 

where q = q{t, x) is a complex-valued function of the two real variables t and x, 
t represents time, and x represents space. q{t, x) is subject to periodic boundary 
condition of period 27r, and even constraint, i.e., 

q{t, X + 2n) = q{t, x) , q{t, —x) — q{t, x) . 

uj e (1/2, 1) is a positive constant, a > and /3 > are constants, and e > is the 
perturbation parameter. The following theorem was proved in |4U) . 

Theorem 2.2. There exists a Eq > 0, such that for any e £ (0,eo); there exists 
a codimension 1 surface in the external parameter space (a, /3, lo) G ]R+ x M+ x 
M+ where lo G (i,l)/S', S is a finite subset, and auj < (3. For any (q;,/3,w) on 
the codimension 1 surface, the singularly perturbed nonlinear Schrddinger equation 
{2.4^ possesses a symmetric pair of Silnikov homoclinic orbits asymptotic to a saddle 
Qe. The codimension 1 surface has the approximate representation given by a = 
1/k{uj), where k{u!) is plotted in Figure^ 
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The proof of the theorem is also completed through two measurements: the Mel- 
nikov measurement and the second measurement. But more powerful machineries 
are needed |4U| . 

2.3. Silnikov Homoclinic Orbits in Vector NLS Under Perturbations. 

In recent years, novel results have been obtained on the solutions of the vector 
nonlinear Schrodinger equations [2] |83| . Abundant ordinary integrable results 
have been carried through 80 21 , including linear stability calculations |20| . 
Specifically, the vector nonlinear Schrodinger equations can be written as 

wt + Pxx + ^{\p\^ + x\q\^)p = 0, 
m + qxx + ^{x\p\^ + \q\'')q^o, 

where p and q are complex valued functions of the two real variables t and x, and 
X is a positive constant. These equations describe the evolution of two orthogonal 
pulse envelopes in birefringent optical fibers |62| j63) . with industrial applications 
in fiber communication systems [25; and all-optical switching devices 28 . For 
linearly birefringent fibers 62 , x = 2/3. For elliptically birefringent fibers, x 
can take other positive values 63 . When % = 1, these equations are first shown 
to be integrable by S. Manakov j5 9j ■ and thus called Manakov equations. When 
X is not 1 or 0, these equations are non-integrable. Propelled by the industrial 
applications, extensive mathematical studies on the vector nonlinear Schrodinger 
equations have been conducted. Like the scalar nonlinear Schrodinger equation, 
the vector nonlinear Schrodinger equations also possess figure eight structures in 
their phase space. Consider the singularly perturbed vector nonlinear Schrodinger 
equations, 

(2.5) ipt +Pxx + \ [{\P?' + kl^) - = ie[Pxx - ap - [3] , 

(2.6) tqt + Qxx + ^[(bP + \q\^) - iO^]q = ~ aq - f3] , 

where p{t, x) and q(t, x) are subject to periodic boundary condition of period 27r, 
and are even in x, i.e. 

p{t, X + 27r) — p{t, x) , p{t, —x) — p{t, x) , 

q(t, x + 2tt) ^ q{t, x) , q{t, ~x) = q{t, x) , 

to e (1, 2), a > and /3 are real constants, and e > is the perturbation parameter. 
We have 

Theorem 2.3 ( '44' ) . There exists a Eq > 0, such that for any e G (0, Eq), there 
exists a codimension 1 surface in the space of {a,f3,uj) G x x where 
U! G (1,2)/S', S is a finite subset, and aoj < \/2/3- -for any {a,(3,U!) on the codi- 
mension 1 surface, the singularly perturbed vector nonlinear Schrodinger equations 
\2. 5|) - i|l^.6)) possesses a homoclinic orbit asymptotic to a saddle Q^. This orbit is 
also the homoclinic orbit for the singularly perturbed scalar nonlinear Schrodinger 
equation studied in last section, and is the only one asymptotic to the saddle Qe for 
the singularly perturbed vector nonlinear Schrodinger equations \2.,^) - \2.f)]) . The 
codimension 1 surface has the approximate representation given by a = l/hi{uj), 
where k{lu) is plotted in Figure^ 



6 



YANGUANG (CHARLES) LI 



Figure 2. The graph of k{uj;N). 



2.4. Silnikov Homoclinic Orbits in Discrete NLS Under Perturba- 
tions . Consider the foUowing perturbed discrete cubic nonhnear Schrodinger equa- 
tions | 5U) . 



iq„ 



(2.7) 

where i = 



1 

+ie 



In+l — 2(7„ + qn-l 



1, g„'s are complex variables, 
Qn+N = Qm (periodic condition); and = g„, (even condition); 



h = jj, and 

TVtan— < w < TVtan— , for iV > 3, 
N N 

TT 

3 tan — < a; < oo, for iV = 3. 
3 

e e [0, £i), a (> 0), /? (> 0) are constants. 

This is a 2{M + 1) dimensional system, where 

(TV even); and M = (iV - l)/2, (TV odd). 

This system is a finite-difference discretization of the perturbed NLS l|2.4|l . The 
following theorem was proved in |50j . 

Denote by Sat (iV > 7) the external parameter space, 

Sat = |(w,a,/3) w e (TVtan ^, iVtan ^), 
a £ (0,ao),/3 e (0,/3o); 

where ao and Pq are any fixed positive numbers. 

Theorem 2.4. For any N (7 < N < oo), there exists a positive number Eq, 
such that for any e € (0, Sq), there exists a codimension 1 surface in T,n; for any 
external parameters (uj,a,P) on E^, there exists a homoclinic orbit asymptotic to a 
saddle Q^. The codimension 1 surface E^ has the approximate expression a = 
where k = k(cl;; N) is shown in Fig\^ 

In the cases (3 < < 6), k is always negative. For N > 7, k can be positive 
as shown in Fig|21 When N is even and > 7, there is in fact a symmetric pair of 
homoclinic orbits asymptotic to a fixed point at the same values of the external 
parameters; since for even N, we have the symmetry: If = f{n,t) solves 1)2. 7|l . 
then qn = f{n + N/2, t) also solves H2.7|l . When N is odd and > 7, the study can 
not guarantee that two homoclinic orbits exist at the same value of the external 
parameters. 
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2.5. Comments on DSII Under Perturbations. Consider the perturbed 
Davey-Stewartson II equations |42j , 



iqt = Tq 



2{\q\'-u') + u, 



q + ief 



Au = -4:dy\q\' 



where g is a complex- valued function of the three variables {t,x,y)^ u is a real- 
valued function of the three variables {t,x,y), T = dxx — dyy, A = dxx + dyy, to > 
is a constant. We also consider periodic boundary conditions. The perturbation / 
can be for example a singular perturbation 

f = Aq-aq + (3 , 

where a > 0, /? > are constants, or a regular perturbation by mollifying A into a 
bounded Fourier multiplier 



Ag = - 2^ /3fe|fc| qk cos kix cos k2y , 
in the case of periods (27r, 27r) 



kez^ 



f3k = 1, <N, (3k ^ \k\-^, \k\ > N, 

for some large N, |fcp = kf + fc|. 

Under both regular and singular perturbations, the rigorous Melnikov mea- 
surement can be established |42j. It turns out that only local well-posedness is 
necessary for rigorously setting up the Melnikov measurement, thanks to the fact 
that the unperturbed homoclinic orbits are classical solutions. The obstacle to- 
ward proving the existence of homoclinic orbits comes from a technical difficulty in 
solving a linear system to get the normal form for proving the size estimate of the 
stable manifold of a saddle. For details, see 1421. 



2.6. Transversal Homoclinic Orbits in a Periodically Perturbed Sine- 
Gordon Equation . Transversal homoclinic orbits in continuous systems often 
appear in two types of systems: (1). periodic systems where the Poincare period 
map has a transversal homoclinic orbit, (2). autonomous systems where the homo- 
clinic orbit is asymptotic to a hyperbolic limit cycle. 

Consider the periodically perturbed sine-Gordon (SG) equation, 

(2.8) utt — (?Uxx + sinu 4- (\—aut + u'^x(||u||) cost], 

where 

1 0, ||u|| > 2M, 

for M < \\u\\ < 2M, x(||u||) is a smooth bump function, under odd periodic bound- 
ary condition, 

u{x + 27r, t) = t), u{x, t) — —u{x, t), 
■| < < 1, a > 0, e is a small perturbation parameter. 

Theorem 2.5 (|lHj, IHSl)- There exists an interval I C M+ such that for any 
a E I, there exists a transversal homoclinic orbit u — S,{x,t) asymptotic to in . 
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2.7. Transversal Homoclinic Orbits in a Derivative NLS . Consider 
the derivative nonlinear Sclirodinger equation, 

(2.9) iqt^qa:x + '2\qfq + ie 

wiiere q is a complex- valued function of two real variables t and x, e > is the 
perturbation parameter, /i is a real constant, and dx is a bounded Fourier multiplier, 

K oo 

Bxq = — kqk sin kx , for q = qk cos kx , 

k=l k=0 

and some fixed K. Periodic boundary condition and even constraint are imposed, 

q{t, X + 2n) ~ q{t, x) , q{t, -x) — q{t, x) . 

Theorem 2.6 (|41j). There exists a Sq > 0, such that for any e E (0, £o); o,nd 
> 5.8, there exist two transversal homoclinic orbits asymptotic to the limit cycle 
Qc = |exp{-i[|t + 7]}. 



3. Existence of Chaos 

The importance of homoclinic orbits with respect to chaotic dynamics was first 
realized by Poincare |68| . In 1961, Smale constructed the well-known horseshoe 
in the neighborhood of a transversal homoclinic orbit |74) |75| |76| . In particular, 
Smale's theorem implies Birkhoff 's theorem on the existence of a sequence of struc- 
turely stable periodic orbits in the neighborhood of a transversal homoclinic orbit 
[7|- In 1984 and 1988 [SS] [521 > Palmer gave a beautiful proof of Smale's theorem 
using a shadowing lemma. Later, this proof was generalized to infinite dimensions 
by Steinlein and Walther |77| (75] and Henry In 1967, Silnikov proved Smale's 
theorem for autonomous systems in finite dimensions using a fixed point argument 
|72) . In 1996, Palmer proved Smale's theorem for autonomous systems in finite 
dimensions using shadowing lemma |15) |K7) . In 2002, Li proved Smale's theorem 
for autonomous systems in infinite dimensions using shadowing lemma |41| . 

For nontransversal homoclinic orbits, the most well-known type which leads 
to the existence of Smale horseshoes is the so-called Silnikov homoclinic orbit |70) 
|71| |73| |17| |18) . Existence of Silnikov homoclinic orbits and new constructions 
of Smale horseshoes for concrete nonlinear wave systems have been established in 
finite dimensions |50| |52j and in infinite dimensions |48| |35| |40| . 

3.1. Shift Automorphism. Let W be a set which consists of elements of the 
doubly infinite sequence form: 

a ~ {■ ■ -0-2(1-1(10, aia2 ■ ■■), 

where Ok G {^i, ■ • • , Im}, ni labels, and k £ Z. We introduce a topology in W by 
taking as neighborhood basis of 

a = {■ ■ ■a_2a_^aQ,a^a2 ■ ■■), 

the set 



Ok = al {\k\ < j) 
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for j = 1,2,- ■ ■. This makes W a topological space. The shift automorphism x is 
defined on W by 

X : W^W, 

Va e W, = ^7 where bk = ak+i- 

The shift automorphism x exhibits sensitive dependence on initial conditions, which 
is a hallmark of chaos. 

3.2. NLS Under Regular Perturbations . We continue from section ITTl 
and consider the regularly perturbed nonlinear Schrodinger (NLS) equation H2.1|l . 
Starting from the Homoclinic Orbit Theorem 12.11 We have the following theorem 
on the existence of chaos. 

Theorem 3.1 (Chaos Theorem, j35j ). Under certain generic assumptions for 
the perturbed nonlinear Schrodinger system 1^2. there exists a compact Cantor 
subset A of (the Sobolev space), A consists of points, and is invariant under a 
Poincare map P. P restricted to A, is topologically conjugate to the shift automor- 
phism X on four symbols 1,2,-1,-2. That is, there exists a homeomorphism 

4> : Wh^A, 

such that the following diagram commutes: 




yv — * A 





Although the symmetric pair of Silnikov homoclinic orbits is not structurally 
stable, the Smale horseshoes are structurally stable. Thus, the Cantor sets and the 
conjugacy to the shift automorphism are structurally stable. 

3.3. NLS Under Singular Perturbations . We continue from section UT^ 
and consider the singularly perturbed nonlinear Schrodinger (NLS) equation ]2A\i . 
Starting from the Homoclinic Orbit Theorem 12. 21 we have 

Theorem 3.2 (Chaos Theorem, 43 ). Under certain generic assumptions for 
the perturbed nonlinear Schrodinger system ^2.4^ , there exists a compact Cantor 
subset A of , A consists of points, and is invariant under a Poincare map P. 
P restricted to A, is topologically conjugate to the shift automorphism x on four 
symbols 1,2,-1,-2. That is, there exists a homeomorphism 

(f> : W^A, 

such that the following diagram commutes: 

W ^ A 

W — > A 
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3.4. Discrete NLS Under Perturbations . We continue from section IT^ 
Starting from tlie Homoclinic Orbit Theorem 12.41 we have the theorem on the 
existence of chaos for the perturbed discrete nonlinear Schrodinger equation H2.7|l . 

Theorem 3.3 (Chaos Theorem, 57). Under certain generic assumptions for 
the perturbed discrete nonlinear Schrodinger system {2.7}) , there exists a compact 
Cantor subset A ofM?^^^^^\ A consists of points, and is invariant under a Poincare 
map P. P restricted to A, is topologically conjugate to the shift automorphism x 
on two symbols 0, 1. That is, there exists a homeomorphism 

(f> : W^A, 

such that the following diagram commutes: 

W ^ A 

X p 

W — * A 



3.5. Numerical Simulation of Chaos. The finite-difference discretization 
of both the regularly and the singularly perturbed nonlinear Schrodinger equations 
(|2.1(l and (|2.4|l leads to the same discrete perturbed nonlinear Schrodinger equation 
lITTIl . 

In the chaotic regime, typical numerical output is shown in Figure 13 Notice 
that there are two typical profiles at a fixed time: (1). a breather type profile with 
its hump located at the center of the spatial period interval, (2). a breather type 
profile with its hump located at the boundary (wing) of the spatial period interval. 
These two types of profiles are half spatial period translate of each other. If we 
label the profiles, with their humps at the center of the spatial period interval, by 
"C"; those profiles, with their humps at the wing of the spatial period interval, by 
"W" , then 

(3.1) 'W"=(To"C", 

where a is the symmetry group element representing half spatial period translate. 
The time series of the output in Figure|31is a chaotic jumping between "C" and "W" , 
which we call "chaotic center-wing jumping" . We interpret the chaotic center-wing 
jumping as the numerical realization of the shift automorphism x on symbols. We 
can make this more precise in terms of the phase space geometry. The figure-eight 
structure of the integrable NLS projected onto the plane of the Fourier component 
cos a: is illustrated in Figure^ and labeled by Lc and Lw- From the symmetry, 
we know that 

Lw = u o Lc- 

Lc has the spatial-temporal profile realization as in Figure|Slwith the hump located 
at the center. L^r corresponds to the half spatial period translate of the spatial- 
temporal profile realization as in Figure El with the hump located at the boundary 
(wing). An orbit inside Lc, Lcin has a spatial-temporal profile realization as in 
Figure ini The half period translate of Lcin, Lwin is inside Lw- An orbit outside 
Lc and Lw, Lout has the spatial-temporal profile realization as in Figure [T] Si 
and Si^a are two phase blocks. From these figures, we can see clearly that the 
chaotic center-wing jumping (Figurc|2Jl is the realization of the shift automorphism 
on symbols in Si U Si^^- 
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Figure 3. A chaotic solution in the discrete perturbed NLS sys- 
tem (j22|l. 

Figure 4. Figure-eight structure and its correspondence with the 
chaotic center- wing jumping. 

Figure 5. Spatial- temporal profile realization of Lc in Figure 01 
(coordinates are the same as in Figure 01 ■ 

Figure 6. Spatial-temporal profile realization of Lcin in Figure 
0] (coordinates are the same as in Figure OJ- 

3.6. Shadowing Lemma and Chaos in Finite-D Periodic Systems . 

Since its invention 3 , shadowing lemma has been a useful tool for solving many 
dynamical system problems. Here we only focus upon its use in proving the exis- 
tence of chaos in a neighborhood of a transversal homoclinic orbit. According to 
the type of the system, the level of difficulty in proving the existence of chaos with 
a shadowing lemma is different. 

A finite-dimensional periodic system can be written in the general form 

X = F{x, t) , 

where x E M", and F{x, t) is periodic in t. Let / be the Poincare period map. 

Definition 3.4. A doubly infinite sequence {j/j} in M" is a (5 pseudo-orbit of 
a diffeomorphism / : R" ^ K" if for all integers j 

\y]+i - fivj)] < <> ■ 

An orbit {f-'{x)} is said to e-shadow the S pseudo-orbit {j/j} if for all integers j 

\rix)-y,\<e. 

Definition 3.5. A compact invariant set S is hyperbolic if there are positive 
constants K, a and a projection matrix valued function P{x), x € S, of constant 
rank such that for all x in S 

P{f{x))Df{x) = Df{x)P{x) , 

\Df^{x)P{x)\<Ke-"\ (j>0), 

\DP{x){I-Pix))\<Ke''^ , (j<0). 

Theorem 3.6 (Shadowing Lemma |66)). Let S be a compact hyperbolic set for 
the diffeomorphism f : i—> R". Then given e > sufficiently small there 
exists S > such that every S pseudo-orbit in S has a unique s -shadowing orbit. 

The proof of this theorem by Palmer |66| is overall a fixed point argument with 
the help of Green functions for linear maps. 

Let yo be a transversal homoclinic point asymptotic to a saddle xq of a 
diffeomorphism / : R" i-^ R". Then the set 

S^{xo}U{r{yo):jeZ} 
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Figure 7. Spatial-temporal profile realization of Lout in Figure 01 
(coordinates are the same as in Figure OJ- 

is hyperbolic. Denote by Aq and Ai the two orbit segments of length 2m + 1 

Ao = {xo,xo,-- - ,xo} , A, {r'"(yo),r'"+'(yo),--- ,r-'(yo),/'"(yo)} • 
Let 

a = (• • • , a_i, ao, ai, • • • ) , 

where aj G {0,1}, be any doubly infinite binary sequence. Let A be the doubly 
infinite sequence of points in S, associated with a 

A = {■ ■ ■ , Aa_i, Aag , Aai , ' ' ' } • 

When m is sufficiently large, A is a 6 pseudo-orbit in S. By the shadowing lemma 
fTheorem l3.fi|l . there is a unique e-shadowing orbit that shadows A. In this manner, 
Palmer |66| gave a beautiful proof of Smale's horseshoe theorem. 

Definition 3.7. Denote by E the set of doubly infinite binary sequences 

«=(■••, a-i, ao, ai, • ■ • ) , 

where aj e {0, 1}. We give the set {0, 1} the discrete topology and S the product 
topology. The Bernoulli shift x is defined by 

[x(a)]j = aj+i . 

Theorem 3.8. Let yo be a transversal homoclinic point asymptotic to a saddle 
Xq of a diffeomorphism f : R" i— > R" . Then there is a homeomorphism <p of T, 
onto a compact subset of R" which is invariant under f and such that when m is 
sufficiently large 

that is, the action of f^™''^^ on is topologically conjugate to the action of x on 
S. 

Here one can define (/>(a) to be the point on the shadowing orbit that shadows 
the midpoint of the orbit segment Aa^, which is cither xq or yo- The topological 
conjugacy can be easily verified. For details, see |66| . Other references can be 
found in |65] |85]. 

3.7. Shadowing Lemma and Chaos in Infinite-D Periodic Systems. 

An infinite-dimensional periodic system defined in a Banach space X can be written 
in the general form 

X = F(x, t) , 

where x G X, and F{x,t) is periodic in t. Let / be the Poincare period map. When 
/ is a map which needs not to be invertible, shadowing lemma and symbolic 
dynamics around a transversal homoclinic orbit can both be established |77| |78| 
j26| . Other references can be found in |24| |13| j86) [8] . There exists also a work on 
horseshoe construction without shadowing lemma for sinusoidally forced vibrations 
of buckled beam |27| . 
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3.8. Periodically Perturbed Sine-Gordon (SG) Equation. We continue 
from section l^in and use the notations in scction rO^ For the periodically perturbed 
sine-Gordon equation H2.8|l . the Poincare period map is a difFeomorphism in H^. 
As a corollary of the result in last section, we have the theorem on the existence of 
chaos. 

Theorem 3.9. There is an integer m and a homeomorphism (j) of T, onto a 
compact Cantor subset A of H^. A is invariant under the Poincare period-2'r: map 
P of the periodically perturbed sine-Gordon equation f2.8\) . The action of P'^'"^^^ 
on A is topologically conjugate to the action of x on T, : p2m+i o cj) = cf) o x- That 
is, the following diagram commutes: 

E ^ A 
E — > A 

3.9. Shadowing Lemma and Chaos in Finite-D Autonomous Systems. 

A finite-dimensional autonomous system can be written in the general form 

i = F{x) , 

where x G K". In this transversal homoclinic orbit can be an orbit asymp- 

totic to a normally hyperbolic limit cycle. That is, it is an orbit in the intersection 
of the stable and unstable manifolds of a normally hyperbolic limit cycle. Instead 
of the Poincare period map as for periodic system, one may want to introduce the 
so-called Poincare return map which is a map induced by the flow on a codimen- 
sion 1 section which is transversal to the limit cycle. Unfortunately, such a map 
is not even well-defined in the neighborhood of the homoclinic orbit. This poses a 
challenging difficulty in extending the arguments as in the case of a Poincare period 
map. In 1996, Palmer |67| completed a proof of a shadowing lemma and existence 
of chaos using Newton's method. It will be difficult to extend this method to infi- 
nite dimensions, since it used heavily differentiations in time. Other references can 
be found in [iSj [121 |Z2] [22] . 

3.10. Shadowing Lemma and Chaos in Inflnite-D Autonomous Sys- 
tems. An infinite-dimensional autonomous system defined in a Banach space X 
can be written in the general form 

i = F{x) , 

where x ^ X. In 2002, the author |41) completed a proof of a shadowing lemma 
and existence of chaos using Fenichel's persistence of normally hyperbolic invariant 
manifold idea. The setup is as follows, 

• Assumption (Al): There exist a hyperbolic limit cycle S and a transver- 
sal homoclinic orbit ^ asymptotic to S. As curves, S and ^ are . 

• Assumption (A2): The Fenichel fiber theorem is valid at S. That is, 
there exist a family of unstable Fenichel fibers {!F^{q) : q £ S} and a 
family of stable Fenichel fibers {J-'^{q) : q £ S}. For each fixed q £ S, 
J^'^iq) and T^q) are submanifolds. and J^%q) are in g,Vg £ 
S. The unions Uges-^"(9) ^^'^ Uqes unstable and stable 
manifolds of S. Both families are invariant, i.e. 

i^*(.F"(g))c.F"(F*(q)),Vi<0,ge5, 
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F'{T'{q))cT%F'{q)),yt>0,qeS, 

where _F* is the evolution operator. There are positive constants k and C 
such that yq G S, Vg" G and yq+ G 

\\F'{q-)-F'{q)\\<Ce^'\\q--q\\yt<0, 

\\F\q+)~F\q)\\<Ce-^'\\q+-q\\,yt>0. 

• Assumption (A3): F*{q) is C° in i, for t G (—00,00), 5 G X. For any 
fixed t G (—00, 00), F*(g) is a diffeomorphism on X. 

Remark 3.10. Notice that we do not assume that as functions of time, S and 
^ are , and we only assume that as curves, S and ^ are C^. 

Under the above setup, a shadowing lemma and existence of chaos can be 
proved |41| . Another crucial element in the argument is the establishment of a 
A-lemma (also called inclination lemma) I41j . 

3.11. A Derivative Nonlinear Schrodinger Equation. We continue from 
section [2.71 and consider the derivative nonlinear Schrodinger equation H2.9|) . The 
transversal homoclinic orbit given in Theorem 12.61 is a classical solution. Thus, 
Assumption (Al) is valid. Assumption (A2) follows from the standard arguments 
in ^52 ^48^ 40 . Since the perturbation in H2.9|l is bounded. Assumption (A3) fol- 
lows from standard arguments. Thus there exists chaos in the derivative nonlinear 
Schrodinger equation H2.9fl (41). 

4. Lax Pairs of Euler Equations of Inviscid Fluids 

The governing equations for the incompressible viscous fluid flow are the Navier- 
Stokes equations. Turbulence occurs in the regime of high Reynolds number. By 
formally setting the Reynolds number equal to infinity, the Navier-Stokes equations 
reduce to the Euler equations of incompressible inviscid fluid flow. One may view 
the Navier-Stokes equations with large Reynolds number as a singular perturbation 
of the Euler equations. 

Results of T. Kato show that 2D Navier-Stokes equations are globally well- 
posed in C°([0,oo);i?^(i?^)), s > 2, and for any < T < 00, the mild solu- 
tions of the 2D Navier-Stokes equations approach those of the 2D Euler equations 
in C°{[0,T];H%R'^)) |3T|. 3D Navier-Stokes equations are locally well-posed in 
C°{[0,t];H''{R^)), s > 5/2, and the mild solutions of the 3D Navier-Stokes equa- 
tions approach those of the 3D Euler equations in C^{[0,t]; H''{R^)), where r de- 
pends on the norms of the initial data and the external force |29| j30j . Extensive 
studies on the inviscid limit have been carried by J. Wu et al. |8H |14j |82j 
There is no doubt that mathematical study on Navier-Stokes (Euler) equations is 
one of the most important mathematical problems. In fact. Clay Mathematics In- 
stitute has posted the global well-posedness of 3D Navier-Stokes equations as one 
of the one million dollars problems. 

V. Arnold |3J realized that 2D Euler equations are a Hamiltonian system. Re- 
cently, the author found Lax pair structures for Euler equations 3_8] |53) |45) |46| 
|51j . Understanding the structures of solutions to Euler equations is of fundamen- 
tal interest. Of particular interest is the question on the global well-posedness of 
3D Navier-Stokes and Euler equations. Our number one hope is that the Lax pair 
structures can be useful in investigating the global well-posedness. Our secondary 
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hope is that the Darboux transformation |53) associated with the Lax pair can 
generate exphcit representation of homocHnic structures |36| . 

The philosophical significance of the existence of Lax pairs for Euler equations is 
even more important. If one defines integrability of an equation by the existence of 
a Lax pair, then both 2D and 3D Euler equations are integrable. More importantly, 
both 2D and 3D Navier-Stokes equations at high Reynolds numbers are singularly 
perturbed integrable systems. Such a point of view changes our old ideology on 
Euler and Navier-Stokes equations. 

4.1. A Lax Pair for 2D Euler Equation. The 2D Euler equation can be 
written in the vorticity form, 

(4.1) dtn + {^,n} = o, 

where the bracket { , } is defined as 

fl is the vorticity, and ^> is the stream function given by, 

u = -dy'i' , V = (9a,* , 
and the relation between vorticity and stream function \1/ is, 

n = dxV- dyu ^ A* . 
Theorem 4.1 (Li, ;38, ). The Lax pair of the 2D Euler equation {4-1}) is given 



Lip = \ip , 



(4.2) 

where 

and X is a complex constant, and ip is a complex-valued function. 

4.2. A Darboux Transformation for 2D Euler Equation. Consider the 
Lax pair H4.2|) at A = 0, i.e. 

(4.3) {n,p} = 0, 

(4.4) dtp + {^,p}^0, 
where we replaced the notation ip by p. 

Theorem 4.2 (!53 ). Let f ~ f(t,x,y) be any fixed solution to the system 
\4-4)) j '^c define the Gauge transform Gf: 

(4.5) p = G/p=^[p.-(a,ln/)p] , 
and the transforms of the potentials 51 and *; 

(4.6) ^ = ^ + F, n = n + AF, 

where F is subject to the constraints 

(4.7) {n, AF} = , {AF, F} = . 

Then p solves the system j^.^ \4-4{ l o.^ (^j^)- Thus \4-5\j and \4-(^ form the 
Darboux transformation for the 2D Euler equation \4-l^ and its Lax pair \4-'A\4-4\ l- 
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Remark 4.3. For KdV equation and many other soliton equations, the Gauge 
transform is of the form |60| . 



In general, Gauge transform does not involve potentials. For 2D Euler equation, a 
potential factor is needed. From H4.3|l . one has 



The Lax pair 14.31 14.4|l has a symmetry, i.e. it is invariant under the transform 
{t,x,y) {—t,y,x). The form of the Gauge transform H4.5I) resulted from the 
inclusion of the potential factor is consistent with this symmetry. 

Our hope is to use the Darboux transformation to generate homoclinic struc- 
tures for 2D Euler equation 36 . 

4.3. A Lax Pair for Rossby Wave Equation. The Rossby wave equation 

is 

dtn + {^i>,n} + pd^'f = , 

where = ^l{t,x,y) is the vorticity, {^t,!]} = "ifx^y ~ '^y^x, and 5* = A~^r2 
is the stream function. Its Lax pair can be obtained by formally conducting the 
transformation, = O + f3y, to the 2D Euler equation 38J, 



where (/? is a complex-valued function, and A is a complex parameter. 

4.4. Lax Pairs for 3D Euler Equation. The 3D Euler equation can be 
written in vorticity form. 



where u — {ui,U2,U3) is the velocity, Q = (fii, f22, i^a) is the vorticity, V = 
{dx-,dy,dz), Q = \7 X u, and V • u = 0. u can be represented by Q for exam- 
ple through Biot-Savart law. 



p = Px- {dx In f)p . 




p ^ _ f^P ^ Py_ _ IilP 



(4.8) 



dtn + {u ■ V)n - (f7 • v)w = , 



Theorem 4.4. The Lax pair of the 3D Euler equation i4-<^ given as 



(4.9) 




where 



L(f)^{n- \/)<f) , A(p={u- V)</> , 
X is a complex constant, and (j) is a complex scalar-valued function. 

Theorem 4.5 ( 10 ). Another Lax pair of the 3D Euler equation j^.^ is given 

as 



(4.10) 




where 



Lip^ {fl- V).^ - {f ■ V)17 , Aip^ {u- V)^ - {if ■ V)u , 



CHAOS IN PDES AND LAX PAIRS OF EULER EQUATIONS 



17 



X is a complex constant, and (p — ((^i, ip2, ^ps) is a complex 3-vector valued Junction. 

Our hope is that the infinitely many conservation laws generated by A G C 
can provide a priori estimates for the global well-posedness of 3D Navier-Stokes 
equations, or better understanding on the global well-posedness |51) . For more 
informations on the topics, sec 53 51 . 

5. Dynamical System Studies on 2D Euler Equation 

5.1. Linearized 2D Euler Equation at a Fixed Point. To begin an in- 
finite dimensional dynamical system study, we investigate the linearized 2D Euler 
equation at a fixed point j37| |45j |45| . We consider the 2D Euler equation (|4.1|l 
under periodic boundary condition in both x and y directions with period 27r. Ex- 
panding Q. into Fourier series, 

O \ " ik-X 

fcez2/{o} 

where w-k — oJk , k — (fci, fc2)^, and X ~ [x, y)^ ■ The 2D Euler equation can be 
rewritten as 

(5.1) iuk^ ^ A{p,q) ujpUjq , 



where A(p, q) is given by, 
(5.2) 



k=p+q 



pi qi 
P2 q2 



where jgp = 91+92 fo^' 9 = (9i; 92)"^, similarly for p. 

Denote {u!k}k£Z^/{o} by w- We consider the simple fixed point u* |37j : 

(5.3) a;; =r, c^* = 0, if fc ^ p or - p, 

of the 2D Euler equation (|5.1|) . where F is an arbitrary complex constant. The 

linearized two-dimensional Euler equation at w* is given by, 

(5.4) ujk = A{p, k-p)T Uk^p + A{-p, k+p)f ujk+p ■ 

Definition 5.1 (Classes). For any k e Z^/{0}, we define the class to be 
the subset of ZV{0}: 



= \ k + npeZy{Q} 



71 e Z, p is specified in (|5.3|) |'. 



See Fig|51for an illustration of the classes. According to the classification defined in 
Definition 15.11 the linearized two-dimensional Euler equation H5.4|l decouples into 
infinitely many invariant subsystems: 

^k+np = Mp. k + in~ l)p) F Wfe+(„_i)p 

(5.5) + A{-^p,k+{n + l)p)f ^i^^^^^^^p. 

Let £^ be the linear operator defined by the right hand side of l|5.5|) . and H'^ be 
the Sobolev space where s > is an integer and 7?" = £2- 
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{-P2,Pl) 




{P2,~Pl) 



Figure 8. An illustration of the classes and the disk D\p\. 



Theorem 5.2 ( |37| 14711. The eigenvalues of the linear operator Cf. in H'^ 
are of four types: real pairs (c,—c), purely imaginary pairs (id, —id), quadruples 
(±c±id), and zero eigenvalues. 

Definition 5.3 (The Disk). The open disk of radius \p\ in Z^/{0}, denoted by 
DIpI, is defined as 



D\p\ = jfc e zV{0} 



\k\ < \P\ 



and the closure of D\p\ , denoted by -D|p| , is defined as 



k e ZV{0} 



\k\ < \P\ 



Theorem 5.4 (The Spectral Theorem, |37] |17|). We have the following claims 
on the spectrum of the linear operator Lj,: 

(1) //E^ n Z)|p| = 0, then the entire spectrum of the linear operator Cf. is 



its continuous spectrum. See Figure\^ where b = — ^|r| 



Pi ki 
P2 k2 



That is, both the residual and the point spectra of Cf, are empty. 
(2) // Ej, n Z?|p| 7^ 0, then the entire essential iJ" spectrum of the linear 
operator Cf. is its continuous spectrum. That is, the residual spectrum of 
Cf. is empty. The point spectrum of Cf is symmetric with respect to both 
real and imaginary axes. See Figure FTTl 

Denote by L the right hand side of (|5.4|) . i.e. the whole linearized 2D Euler operator, 
the spectral mapping theorem holds. 
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3{A} 





i2\b\ 




* 




-i2\b\ 



5R{A} 



Figure 9. The spectrum of £^ in case (1). 



Theorem 5.5 ( 34 ). 

Let C denote the number of points q G Z^/{0} that belong to the open disk of 
radius \p\, -D|p|, and such that q is not paraUel to p. 

Theorem 5.6 ( 34 ). The number of nonimaginary eigenvalues of L (counting 
the multiplicities) does not exceed 2^. 

Another interesting discussion upon the discrete spectrum can be found in |19| . 

Since the introduction of continued fractions for calculating the eigenvalues 
of steady fluid flow, by Meshalkin and Sinai |64) . this topics had been extensively 
explored |H1| ED EH HZ] HH] HI I3Z1- Rigorousjustification on the continued 
fraction calculation was given in |37| |58| . As an example, we take p = (1,1)-^. 
When r ^ 0, the fixed point has 4 eigenvalues which form a quadruple. These four 
eigenvalues appear in the invariant linear subsystem labeled by fc = (—3, —2)^. 
One of them is j371 : 

(5.6) A = 2A/|r| = 0.24822302478255 + i 0.35172076526520 . 

See Figure [TTl for an illustration. The essential spectrum (= continuous spectrum) 
of Cj. with k = (—3, —2)-^ is the segment on the imaginary axis shown in Figure [TTl 
where b = The essential spectrum {— continuous spectrum) of the linear 2D 

Euler operator at this fixed point is the entire imaginary axis. A rather well-known 
open problem is proving the existence of unstable, stable, and center manifolds. The 
main difficulty comes from the fact that the nonlinear term is non-Lipschitzian. 
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i2\b\ 


• 


• 


• 


• 


• 


• 




* 


• 


• 


• 


• 


• 






-i2\b\ 


Figure 10. The spect 


rum of Cj^ in case (2). 



5R{A} 



5.2. Models. To simplify our study, we study only the case when ujk is real, 
Vfc e 1? /{O}, i.e. we only study the cosine transform of the vorticity, 

il = ^ uJk cos(fc • X) . 

fc«EZ2/{0} 

To further simplify our study, we will study a concrete dashed-line model based 
upon the line of fixed points (|5.3|l with the mode p = (1, 1)"^ parametrized by F, 
and the only unstable invariant linear subsystem labeled by fc = (~3, —2)'^. The 
so-called dashed-line model is given by [39) . 

(5.7) 
where 

(5.8) ujn = ^k+np ' ^« = ^(P' k + np) , Am,n = + mp, k + np) , 

^ r 1 , if n ^ 5j ,yjez , 

" \ e , if n = 5j , for some j E Z . 

The model is designed to model the hyperbolic structure of the 2D Euler equation, 
connected to the line of fixed points H5.3|l with p = (1, 1)-^. FigurcIT^ illustrates the 
collocation of the modes in this model, which has the "dashed-line" nature leading 
to the name of the model. When e = 0, one can determine the stable and unstable 
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U{X} 



-i2\b\ 



Figure 11. The spectrum of £f, with k = (-3, -2)"^, when p = (1, 1)^. 
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Figure 12. The collocation of the modes in the dashed- line model, 
manifolds of the fixed point iv*: 

(5.9) ujp = T,(Jn = (n e Z) . 
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Figure 13. The unstable and stable manifolds of the fixed point 
uj* (|5.3(l : (a) and (b) show two "painted eggs" and (c) shows a 
"hp" . 



Using the polar coordinates: 

uji = r cos , UJ4 = r sin 6 ; UJ2 = p cos •& , 103 = p sin ■& ; 

we have the following explicit expressions for the stable and unstable manifolds of 
the fixed point to* and its negative —uj* represented through perverted heteroclinic 
orbits: 



r tanh T , 



Ai 



r sech T , 



(5.10) 



— - In cosh T 
2k 



A. 



arcsm 



?9 = 



A2 



arcsm 



A2 

-Ai 



{k > 0) 
(k < 0) 



where r = nVt+To, (tq, Oq) are the two parameters parametrizing the two-dimensional 
stable (unstable) manifold, and 



K = V-^1^2 COs{e + l9) = ±^/~AlA2^1 

The two auxilliary variables ujo and loq have the expressions: 



A2_ 

AAi 



UJQ 

where 



af3 



sech T ■|^sin[/3 In cosh r - 
sech r cos[/3 In coshr 



a — —AiTk, 



Ai 



cos[/3 In coshr 
sin[/3 In coshr 

2k 



A2 - Ai 

The explicit expression (|5.1U|) of the perverted heteroclinic orbits shows that the 
unstable manifold H^"(w*) of the fixed point w* is the same as the stable manifold 
W^{ijj*) of the negative —lo* of uj* , and the stable manifold of the fixed point lo* 
is the same as the unstable manifold of the negative —lo* of lo* . Both W^{u)*) 
and ^¥"[1^*) have the shapes of "painted eggs" (Figure [T^ . VF"(aj*) and W^{lo*) 
together form the "lip" (Figure which is a higher dimensional generalization of 
the heteroclinic connection on plane. 

As the first step toward understanding the degeneracy v.s. nondegeneracy of 
the hyperbolic structure of 2D Euler equation, we are interested in the e-homotopy 
deformation of such hyperbolic structure for the dashed-line model. We have cal- 
culated the Melnikov functions to study the "breaking" or "persistence" of such 
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hyperbolic structure when e is small. It turns out that both the first and the 
second order Melnikov functions are identically zero |39| . 

6. Conclusion 

Results on chaos in partial differential equations are reported. Results on Lax 
pairs of Euler equations of incompressible inviscid fluids are also reported. Prelim- 
inary results on dynamical system studies of 2D Euler equation are summarized. 
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